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Abstract. We describe new families of random fractals, referred to as "In- 
variable" , which are intermediate between the notions of deterministic and of 
standard random fractals. The parameter V describes the degree of "variabil- 
ity": at each magnification level any variable fractals has at most V key 
"forms" or "shapes" . V-variable random fractals have the surprising property 
that they can be computed using a forward process. More precisely, a version of 
the usual Random Iteration Algorithm, operating on sets (or measures) rather 
than points, can be used to sample each family. To present this theory, we 
review relevant results on fractals (and fractal measures), both deterministic 
and random. Then our new results are obtained by constructing an iterated 
function system (a super IFS) from a collection of standard IFSs together with 
a corresponding set of probabilities. The attractor of the super IFS is called a 
superfractal; it is a collection of V-variable random fractals (sets or measures) 
together with an associated probability distribution on this collection. When 
the underlying space is for example K^, and the transformations are compu- 
tationally straightforward (such as affine transformations), the superfractal 
can be sampled by means of the algorithm, which is highly efficient in terms 
of memory usage. The algorithm is illustrated by some computed examples. 
Some variants, special cases, generalizations of the framework, and potential 
applications are mentioned. 



1. Introduction and Notation 

1.1. Fractals and Random Fractals. A theory of deterministic fractal sets and 
measures, using a "backward" algorithm, was developed in Hutchinson |1()| . A 
different approach using a "forward" algorithm was developed in Barnsley and 
Demko g]. 

Falconer jjl], Graf jjl] and Mauldin and Williams [21] randomized each step in 
the backward construction algorithm to obtain random fractal sets. Arbeiter [Jj 
introduced and studied random fractal measures; see also Olsen [33] ■ Hutchinson 
and Riichendorff jl7| and |18) introduced new probabilistic techniques which allowed 
one to consider more general classes of random fractals. For further material see 
Zahle [5S], Patzschke and Zahle and the references in all of these. 

This paper begins with a review of material on deterministic and random fractals 
generated by IFSs, and then introduces the class of F-variable fractals which in a 
sense provides a link between deterministic and "standard" random fractals. 
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Deterministic fractal sets and measures are defined as tfie attractors of certain 
iterated function systems (IFSs), as reviewed in Section|21 Approximations in prac- 
tical situations quite easily can be computed using the associated random iteration 
algorithm. Random fractals are typically harder to compute because one has to 
first calculate lots of fine random detail at low levels, then one level at a time, build 
up the higher levels. 

In this paper we restrict the class of random fractals to ones that we call random 
^-variable fractals. Superfractals are sets of F-variable fractals. They can be 
defined using a new type of IFS, in fact a "super" IFS made of a finite number N 
of IFSs, and there is available a novel random iteration algorithm: each iteration 
produces new sets, lying increasingly close to V-variable fractals belonging to the 
superfractal, and moving ergodically around the superfractal. 

Superfractals appear to be a new class of geometrical object, their elements ly- 
ing somewhere between fractals generated by IFSs with finitely many maps, which 
correspond to V = N = 1, and realizations of the most generic class of random 
fractals, where the local structure around each of two distinct points are indepen- 
dent, corresponding to y = oo. They seem to allow geometric modelling of some 
natural objects, examples including realistic-looking leaves, clouds, and textures; 
and good approximations can be computed fast in elementary computer graphics 
examples. They are fascinating to watch, one after another, on a computer screen, 
diverse, yet ordered enough to suggest coherent natural phenomena and potential 
applications. 

Areas of potential applications include computer graphics and rapid simulation 
of trajectories of stochastic processes The forward algorithm also enables rapid 
computation of good approximations to random (including "fully" random) pro- 
cesses, where previously there was no available efficient algorithm. 

1.2. An Example. Here we give an illustration of an application of the theory in 
this paper. By means of this example we introduce informally V-variable fractals 
and superfractals. We also explain why we think these objects are of special interest 
and deserve attention. 

We start with two pairs of contractive affine transformations, {/i,/2} and 
where : □ ^ □ with □ := [0,1] x [0,1] C We use two pairs 
of screens, where each screen corresponds to a copy of □ and represents for exam- 
ple a computer monitor. We designate one pair of screens to be the Input Screens, 
denoted by (□i,n2). The other pair of screens is designated to be the Output 
Screens, denoted by (□i',n2')- 

Initialize by placing an image on each of the Input Screens, as illustrated in 
Figure 12 and clearing both of the Output Screens. We construct an image on each 
of the two Output Screens as follows. 

(i) Pick randomly one of the pairs of functions {/i , /I } or {/i j /I }j say{/"\ /^^l- 
Apply f^^ to one of the images on Di or selected randomly, to make an image 
on Di'. Then apply /^^ to one of the images on Di or ^2, also selected randomly, 
and overlay the resulting image / on the image now already on Di/. (For example, 
if black-and-white images are used, simply take the union of the black region of / 
with the black region on Di/, and put the result back onto Di'.) 

(ii) Again pick randomly one of the pairs of functions {/i , /2} or {/i , /I }, say 
{/r^'/?^}- Apply fi^ to one of the images on Di, or ^2, selected randomly, to 
make an image on 02'. Also apply /|^^ to one of the images on Di, or ^2, also 
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Figure 1. Triangles used to define the four transformations 
/i,/i,/2,and/|. 



selected randomly, and overlay the resulting image on the image now already on 
□2'. 

(iii) Switch Input and Output, clear the new Output Screens, and repeat steps 
(i), and (ii). 

(iv) Repeat step (iii) many times, to allow the system to settle into its "stationary 
state" . 

What kinds of images do we see on the successive pairs of screens, and what are 
they like in the "stationary state" ? What does the theory developed in this paper 
tell us about such situations? 

As a specific example, let us choose 



(1.1) 
(1.2) 



1, > ,1 3 51 3 3> 
/i(-,y) = (2^-8y+Y6'r+8^+16)' 



^1, , ,1 3 3 13 11, 

/2(-,y) = (2-+8j/ + Y^,-2-+8y + Y^), 

^2, ^ ,1 3 5 13 13, 
fi (x.v) — —X ?/ H , X y H , 



(1.3) 



3 1 3 5 , 
tx- -y+ —)■ 



1 

2^" 8" ' 16' 2" 8" ' 16' 
We describe how these transformations act on the triangle ABC in the diamond 
ABCD, where A = (i, i), B = (i, |), C = (f , i), and D = (i, i). Let B, = 
(i' D' ^2 = (f, f), B3 - (^, ^), and S4 = (f , f ). See Figureffl Then we 
have 

fl{A) = A, fl{B) ^ B^, fliC) = B; 
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Figure 3. The pair of images after one iteration. 

f^{A)=B,f^iB)=B2,n{C)^C- 

fl{A) = A, fl{B) ^ B^, fl{C) ^ D- 

fl{A)^D,fiiB)=Bi,fiiC)^C. 
In Figure 121 we show an initial pair of images, two jumping fish, one on each 
of the two screens Di and In Figures |31 EI [3 [HI [3 |H1 and we show the 
start of the sequence of pairs of images obtained in a particular trial, for the first 
seven iterations. Then in Figures E| ^2 ^^nd ^1 we show three successive pairs 
of computed screens, obtained after more than twenty iterations. These latter 
images are typical of those obtained after twenty or more iterations, very diverse, 
but always representing continuous "random" paths in R^; they correspond to the 
"stationary state" , at the resolution of the images. More precisely, with probability 
one the empirically obtained distribution on such images over a long experimental 
run corresponds to the stationary state distribution. 

Notice how the two images in Figure lTTl consist of the union of shrunken copies of 
the images in Figure ITUl while the curves in Figure [T^ are made from two shrunken 
copies of the curves in FigurelTTl 
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Figure 6. The two images after four iterations. Both images are 
again the same, a braid of fish. 
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Figure 7. The two images after five iterations. The two images 
are the same. 




Figure 9. The two images after seven iterations. 
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Figure 12. Images on the two screens after L + 2 iterations. 
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This example illustrates some typical features of the theory in this paper, (i) New 
images are generated, one per iteration per screen, (ii) After sufficient iterations 
for the system to have settled into its "stationary state", each image looks like a 
finite resolution rendering of a fractal set that typically changes from one iteration 
to the next; each fractal belongs to the same family, in the present case a family 
of continuous curves, (iii) In fact, it follows from the theory that the pictures in 
this example correspond to curves with this property: for any e > the curve is 
the union of "little" curves, ones such that the distance apart of any two points 
is no more than e, each of which is an affine transformation of one of at most 
two continuous closed paths in R^. (iv) We will show that the successive images, 
or rather the abstract objects they represent, eventually all lie arbitrarily close to 
an object called a superfractal. The superfractal is the attractor of a superlFS 
which induces a natural invariant probability measure on the superfractal. The 
images produced by the algorithm are distributed according to this measure, (v) 
The images produced in the "stationary state" are independent of the starting 
images. For example, if the initial images in the example had been of a dot or a 
line instead of a fish, and the same sequence of random choices had been made, 
then the images produced in Figures EH and El would have been the same at 
the printed resolution. 

One similarly obtains V- variable fractals and their properties using V, rather 
than two, screens and otherwise proceeding similarly. In (iii) each of the sets of 
diameter at most e is an afhne transformation of at most V sets in M^, where these 
sets again depend upon e and the particular image. 

This example and the features just mentioned suggest that superfractals are of 
interest because they provide a natural mathematical bridge between deterministic 
and random fractals, and because they may lead to practical applications in digital 
imaging, special effects, computer graphics, as well as in the many other areas where 
fractal geometric modelling is applied. 

1.3. The structure of this paper. The main contents of this paper, while con- 
ceptually not very difficult, involves potentially elaborate notation because we deal 
with iterated function systems (IFSs) made of IFSs, and probability measures on 
spaces of probability measures. So a material part of our effort has been towards 
a simplified notation. Thus, below, we set out some notation and conventions that 
we use throughout. 

The core machinery that we use is basic IFS theory, as described in and 
0. So in Section El we review relevant parts of this theory, using notation and 
organization that extends to and simplifies later material. To keep the structural 
ideas clear, we restrict attention to IFSs with strictly contractive transformations 
and constant probabilities. Of particular relevance to this paper, we explain what 
is meant by the random iteration algorithm. We illustrate the theorems with simple 
applications to two-dimensional computer graphics, both to help with understand- 
ing and to draw attention to some issues related to discretization that apply a 
fortiori in computations of V-variable fractals. 

We begin Section with the definition of a superlFS, namely an IFS made of 
IFSs. We then introduce associated trees, in particular labelled trees, the space 
of code trees J7, and construction trees; then we review standard random fractals 
using the terminology of trees and superlFSs. 
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In Section 0] we study a special class of code trees, called F-variable trees, 
where V is an integer. What are these trees like? At each level they have at 
most V distinct subtrees! In fact these trees are described with the aid of an IFS 
{n^;7f,r'',ae A} where ^ is a finite index set, V^s are probabilities, and each ry" 
is a contraction mapping from fl^ to itself. The IFS enables one to put a measure 
attractor on the set of ^-variable trees, such that they can be sampled by means 
of the random iteration algorithm. We describe the mappings 77" and composi- 
tions of them using certain finite doubly labelled trees. This, in turn, enables us 
to establish the convergence, as ^ ^ cxd, of the probability measure on the set of 
^-variable trees, associated with the IFS and the random iteration algorithm, to a 
corresponding natural probability distribution on the space fl. 

In Section 13 the discussion of trees in Section 01 is recapitulated twice over: the 
same basic IFS theory is applied in two successively more elaborate settings, yielding 
the formal concepts of V-variable fractals and superfractals. More specifically, in 
Sectional, the superlFS is used to define an IFS of functions that map I/-tuples of 
compact sets into V-tuples of compact sets; the attractor of this IFS is a set of V- 
tuples of compact sets; these compact sets are named T^-variable fractals and the set 
of these variable fractals is named a superfractal. We show that these T^- variable 
fractals can be sampled by means of the random iteration algorithm, adapted to the 
present setting; that they are distributed according to a certain stationary measure 
on the superfractal; and that this measure converges to a corresponding measure 
on the set of "fully" random fractals as V 00, in an appropriate metric. We 
also provide a continuous mapping from the set of ^-variable trees to the set of 
^-variable fractals, and characterize the ^-variable fractals in terms of a property 
that we name "T^- variability" . Section [S] 2 follows the same lines as in Section 
HJl, except that here the superlFS is used to define an IFS that maps l^-tuples 
of measures to V^-tuples of measures; this leads to the definition and properties of 
T^-variable fractal measures. In Sectionals we describe how to compute the fractal 
dimensions of V-variable fractals in certain cases and compare them, in a case 
involving Sierpinski triangles, with the fractal dimensions of deterministic fractals, 
"fully" random fractals, and "homogeneous" random fractals that correspond to 
V — I and are a special case of a type of random fractal investigated by Hambly 
and others 151, 0, [El, 123 • 

In Sectioniniwe describe some potential applications of the theory including new 
types of space-filling curves for digital imaging, geometric modelling and texture 
rendering in digital content creation, and random fractal interpolation for computer 
aided design systems. In SectionjTlwe discuss generalizations and extensions of the 
theory, areas of ongoing research, and connections to the work of others. 

1.4. Some Notation. We use notation and terminology consistent with 

Throughout we reserve the symbols M, N, and V for positive integers. We will 

use the variables m e {1,2, ...,M}, n G {1,2, ...,N}, and v e {1,2, V}. 

Throughout we use an underlying metric space (X, dx) which is assumed to be 

compact unless otherwise stated. We write to denote the compact metric space 

X X X X ... X X. 

V TIMES 

with metric 

d{x,y) = dxv{x,y) = max {dxixy,yv) | w = 1,2, ...,V}, \fx,y G X^, 
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where x = {xi,X2, ■■■,xv) and y = {yi,y2, ■■■,yv)- 

In some applications, to computer graphics for example, (X, o?x) is a bounded 
region in with the Euclidean metric, in which case we will usually be concerned 
with affine or projective maps. 

Let § = S(X) denote the set of all subsets of X, and let C G §. We extend the 
definition of a function / :X^Xto/:S^Sby 

/(C) = {fix) \xeC} 

Let H=H(X) denote the set of non-empty compact subsets of X. Then if / : X — > 
X we have / : HI — > H. We use dm to denote the Hausdorff metric on H implied 
by the metric dx on X. This is defined as follows. Let A and B be two sets in H, 
define the distance from A to B to be 

(1.4) V{A, B) = max{min{dx(a;, y) \ y ^ B} \ x A], 

and define the Hausdorff metric by 

dn{A, B) = max{X'(74, B),V{B, A)}. 

Then (IHI,dH) is a compact metric space. We will write (H^, dgv) to denote the V- 
dimensional product space constructed from (IHI,(iH) just as (X^, dxv ) is constructed 
from (X, dx)- When we refer to continuous, Lipschitz, or strictly contractive func- 
tions acting on we assume that the underlying metric is d^v . 

We will in a number of places start from a function acting on a space, and extend 
its definition to make it act on other spaces, while leaving the symbol unchanged 
as above. 

Let B = B(X) denote the set of Borel subsets of X. Let P = P(X). In some 
applications to computer imaging one sets X = [0, 1] x [0, 1] C and identifies 
a black and white image with a member of EII(X). Greyscale images are identified 
with members of P(X). Probability measures on images are identified with P(]HI(X)) 
or P(P(X)). 

Let (ip(x) denote the Monge Kantorovitch metric on P(X). This is defined as 
follows. Let /i and ly be any pair of measures in P. Then 

dp(/i,i') = sup fdfi~J fdiy\f:X^R, \f{x)-f{y)\<dx{x,y)yx,ye. 

Then (P, dp) is a compact metric space. The distance function dp metrizes the 
topology of weak convergence of probability measures on X, We define the 
push-forward map / : P(X) ^ P(X) by 



Again here we have extended the domain of action of the function / : X ^ X. 

We will use such spaces as P(H^) and P((P(X))^) (or H(H^) and H(P^) de- 
pending on the context). These spaces may at first seem somewhat Baroque, but 
as we shall see, they are very natural. In each case we assume that the metric of 
a space is deduced from the space from which it is built, as above, down to the 
metric on the lowest space X, and often we drop the subscript on the metric without 
ambiguity. So for example, we will write 

d{A, B) = dH((P(x))-) (A B) V A S e H((P(X))^). 



We also use the following common symbols: 
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N= {1,2,3,...}, Z = {... -2,-1,0,1,2,...}, and Z+= {0,1,2,...}. 
When S* is a set, IS*] denotes the number of elements of S. 

2. Iterated Function Systems 

2.1. Definitions and Basic Results. In this section we review relevant infor- 
mation about IFSs. To clarify the essential ideas we consider the case where all 
mappings are contractive, but indicate in Sectional how these ideas can be gener- 
alized. The machinery and ideas introduced here are applied repeatedly later on in 
more elaborate settings. 
Let 

(2.1) {X; fi, f2, ■■; fM;Pi,P2, ■■■,Pm} 

denote an IFS with probabilities. The functions /,„ : X — > X are contraction map- 
pings with fixed Lipschitz constant < Z < 1; that is 

difmix). fra{y)) < I ' d{x, y) Vx, y G X, Vm G {1, 2, M}. 

The Pm's are probabilities, with 

M 

^ = 1, p,n > Vm. 

m—1 

We define mappings F : M(X) H(X) and F : P(X) ^ P(X) by 

M 

F{K) = y UXK) yK e H, 

m— 1 

and 

M 

F{fi) - PmfM V/i e p. 

m—1 

In the latter case note that the weighted sum of probability measures is again a 
probability measure. 

Theorem 1. [HI The mappings F : H(X) H(X) and F : P(X) P(X) are both 
contractions with factor < I < 1. That is, 

d{F{K), F{L)) < I ■ d{K, L) V K,Le H(X), 

and 

d{F{fi),F{,y)) < I ■ d{n, P(X). 
As a consequence, there exists a unique nonempty compact set A G IHI(X) such that 

F{A) = A, 
and a unique measure fi G P(X) such that 

F{p) - 

The support of fi is contained in, or equal to A, with equality when all of the 
probabilities Pm are strictly positive. 

Definition 1. The set A in Theorem^is called the set attractor of the IFS F, and 

the measure fx is called the measure attractor of F. 
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We will use the term attractor of an IFS to mean either the set attractor or the 
measure attractor. We will also refer informally to the set attractor of an IFS as a 
fractal and to its measure attractor as a fractal measure, and to either as a fractal. 
Furthermore, we say that the set attractor of an IFS is a deterministic fractal. This 
is in distinction to random fractals, and in particular to V -variable random fractals 
which are the main goal of this paper. 

There are two main types of algorithms for the practical computation of attrac- 
tors of IFS that we term deterministic algorithms and random iteration algorithms, 
also known as backward and forward algorithms, c.f. (8j. These terms should not 
be confused with the type of fractal that is computed by means of the algorithm. 
Both deterministic and random iteration algorithms may be used to compute deter- 
ministic fractals, and as we discuss later, a similar remark applies to our ^-variable 
fractals. 

Deterministic algorithms are based on the following: 
Corollary 1. Let Aq G H(X), or /Iq G P(X), and define recursively 
Ak = F{Ak-i), or fi,, = Fifif,_^), Vfc G N, 

respectively; then 

(2.2) lim Ak = A, or lim fi^. = /i, 

k — ^oo k — >oo 

respectively. The rate of convergence is geometrical; for example, 

d{Ak,A) < I'' ■d{Ao,A) Vfc G N. 

In practical applications to two-dimensional computer graphics, the transfor- 
mations and the spaces upon which they act must be discretized. The precise 
behaviour of computed sequences of approximations to an attractor of an IFS de- 
pends on the details of the implementation and is generally quite complicated; for 
example, the discrete IFS may have multiple attractors, see Chapter 4. The 
following example gives the flavour of such applications. 

Example 1. In Figure [T^ we illustrate a practical deterministic algorithm, based 
on the first formula in Equation (|2.2|l starting from a simple IFS on the unit square 
□ C R^. The IFS is = {□; fl, fl, /|; 0.36, 0.28, 0.36} where the transformations 
are defined in Equations Hl.l|) . H1.2|l . and (|1.3|) . The successive images, from left 
to right, from top to bottom, represent Aq, A2, A<^, Aj, A20, and A21. In the 
last two images, the sequence appears to have converged at the printed resolution 
to representations of the set attractor. Note however that the initial image is 
partitioned into two subsets corresponding to the colours red and green. Each 
successive computed image is made of pixels belonging to a discrete model for □ 
and consists of red pixels and green pixels. Each pixel corresponds to a set of points 
in ffi.^. But for the purposes of computation only one point corresponding to each 
pixel is used. When both a point in a red pixel and point in a green pixel belonging 
to say An are mapped under F to points in the same pixel in An+i a choice has to 
be made about which colour, red or green, to make the new pixel of An+i. Here 
we have chosen to make the new pixel of A„^i the same colour as that of the pixel 
containing the last point in A„, encountered in the course of running the computer 
program, to be mapped to the new pixel. The result is that, although the sequence 
of pictures converge to the set attractor of the IFS, the colours themselves do not 
settle down, as illustrated in Figure El We call this "the texture efi^ect" , and 
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Figure 13. An illustration of the deterministic algorithm. 

comment on it in Example 13 In printed versions of the figures representing A20 , 
and A21 the red and green pixels are somewhat blended. 

The following theorem is the mathematical justification and description of the 
random iteration algorithm. It follows from Birkhoff's ergodic theorem and our 
assumption of contractive maps. A more general version of it is proved in |l(Jj . 

Theorem 2. Specify a starting point xi e X. Define a random orbit of the IFS to 
be {xi}f^^ where xi+i — fm{xi) with probability Pm ■ Then for almost all random 
orbits {xi}'j^^ we have: 

(2.3) in^K^^i^il^illZlM. 

/ — >oo I 

for all B G B(X) such that fJ-{dB) — 0, where dB denotes the boundary of B. 

Remark 1. This is equivalent by standard arguments to the following: for any 
xi G X and almost all random orbits the sequence of point measures j{Sxi + 5x2 + 
... + 6x1) converges in the weak sense to pi, see for example [?], pages 11 and 12. 
(Weak convergence of probability measures is the same as convergence in the Monge 
Kantorovitch metric, see 0, pages 310 and 311.) 

The random iteration algorithm can be applied to the computation of two- 
dimensional computer graphics. It has benefits compared to deterministic iteration 
of low memory requirements, high accuracy — as the iterated point can be kept at 
much higher precision than the resolution of the computed image — and it allows 
the efficient computation of zooms into small parts of an image. However, as in 
the case of deterministic algorithms, the images produced depend on the computa- 
tional details of image resolution, the precision to which the points {xi,X2^ 
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Figure 14. "Picture" of the measure attractor of an IFS with 
probabilities produced by the random iteration algorithm. The 
measure is depicted in shades of green, from (black) to 255 (bright 
green) . 



are computed, the contractivity of the transformations, the way in which Equation 
(12.311 is implemented, choices of colours, etc. Different implementations can produce 
different results. 

Example 2. Figure El shows a "picture" of the invariant measure of the IFS 
in Example n computed using a discrete implementation of the random iteration 
algorithm, as follows. Pixels corresponding to a discrete model for □ C are 
assigned the colour white. Successive floating point coordinates of points in □ are 
computed by random iteration and the first (say) one hundred points are discarded. 
Thereafter, as each new point is calculated, the pixel to which it belongs is set to 
black. This phase of the computation continues until the pixels cease to change, 
and produces a black image of the support of the measure, the set attractor of the 
IFS, against a white background. Then the random iteration process is continued, 
and as each new point is computed the green component of the pixel to which the 
latest point belongs is brightened by a fixed amount. Once a pixel is at brightest 
green, its value is not changed when later points belong to it. The computation 
is continued until a balance is obtained between that part of the image which is 
brightest green and that which is lightest green, and then stopped. 

The following theorem expresses the ergodicity of the IFS F. The proof depends 
centrally on the uniqueness of the measure attractor. A variant of this theorem, 
weaker in the constraints on the IFS but stronger in the conditions on the set B, 
and stated in the language of stochastic processes, is given by j^. We prefer the 
present version for its simple statement and direct measure theoretic proof. 

Theorem 3. Suppose that fi is the unique measure attractor for the IFS F. Suppose 
B e B(X) is such that f,n{B) C B ym e {1, 2, A/}. Then n{B) = or 1. 
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Proof. Let us define the measure i.i[B (/i restricted by B) by {fj.[B){E) — ^{Br\E). 
The main point of the proof is to show that /i[-B is invariant under the IFS F. ( A 
similar result applies to where B'-^ denotes the complement of B.) 

li E d B^ , for any m, since fm{B) C B, 

fMB){E) = fi{B n f,-\E)) = /i(0) = 0. 

Moreover, 

(2.4) ti{B) = fMBm) - /„(/iLi?)(i3). 
It follows that 

M M M 

fi{B) = PmfmlJi{B) = ^ P™/„(/4S)(S) + E Pm/m(A^LS'')(^) 

m—1 m—1 m—1 

A/ 

= + ^ Pm/m(AiL^^)(S) (from dH). 

m=l 

Hence 

M 

(2.5) Y.P^r.fMB^){B)=0. 

m—1 

Hence for any measurable set C X 

M 

(AiLs)(s) = fi{B n E) = PmfmKB n E) 

m—1 

M M 



= J2 PmfMB){B r)E)+Y, PrnfMB^){B H E) 
m—1 m—1 
A/ 

= J2 PmfMBKE) + (using (ESJ). 



Thus ^ [-B is invariant and so is either the zero measure or for some constant c > 1 
we have cii[B = fi (by uniqueness) = fi[B + ix\_B^ . This implies /i[-B'^ — and in 
particular ii{B'-^) — and ^l{B) = 1. | 

Example 3. Figure El shows close-ups on the two images at the bottom left in 
Figure E| see Example ^ At each iteration it is observed that the pattern of red 
and green pixels changes in a seemingly random manner. A similar texture effect is 
often observed in other implementations and in applications of y- variable fractals to 
computer graphics. Theorem 13 provides a simple model explanation for this effect 
as follows. Assume that the red pixels and the green pixels both correspond to sets 
of points of positive measure, both invariant under F. Then we have a contradiction 
to the Corollary above. So neither the red nor the green set can be invariant under 
F. Hence, either one of the sets disappears — which occurs in some other examples 
— or the pixels must jump around. A similar argument applied to powers of F 
shows that the way the pixels jump around cannot be periodic, and hence must be 
"random" . (A more careful explanation involves numerical and statistical analysis 
of the specific computation.) 
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Figure 15. Close-up on the same small region in each of the bot- 
tom left two images in Figure El showing the texture effect; the 
distribution of red and green pixels changes with each iteration. 

2.2. Fractal Dimensions. In the literature there are many different definitions of 
a theoretical quantity called the "fractal dimension" of a subset of X. A mathemat- 
ically convenient definition of the fractal dimension of a set S* C X is the Hausdorff 
dimension. This is always well-defined. Its numerical value often but not always 
coincides with the values provided by other definitions, when they apply. 

Fractal dimensions are useful for a number of reasons. They can be used to 
compare and classify sets and measures and they have some natural invariance 
properties. For example the Hausdorff dimension of a set 5* is invariant under any 
bi-Lipshitz transformation; that is, if / : X — > X is such that there are constants 
Ciand C2 in (0, oo) with ci • d(x,y) < d{f{x),f{y)) < c\ ■ d{x,y) V x,?/ G X then 
the Hausdorff dimension of S is the same as that of f{S). Fractal dimensions are 
useful in fractal image modelling: for example, empirical fractal dimensions of the 
boundaries of clouds can be used as constraints in computer graphics programs for 
simulating clouds. Also, as we will see below, the specific value of the Hausdorff 
dimension of the set attractor A of an IFS can yield the probabilities for most 
efficient computation of A using the random iteration algorithm. For these same 
reasons fractal dimensions are an important concept for V-variable fractals and 
superfractals. 

The following two definitions are discussed in [l^ pp. 25 et seq. 
Definition 2. Let S cX, S > 0, and < s < oo. Let 



where If/ij" denotes the s"* power of the diameter of the set Ui, and where a 5— cover 
of S is a covering of S by subsets of X of diameter less than d. Then the s- 
dimensional Hausdorff measure of the set S is defined to be 




H'iS) = limm 
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The s-dimensional Hausdorff measure is a Borel measure but is not normally 
even cr-finite. 

Definition 3. The Hausdorff dimension of the set 5 C X is defined to be 

diniH S = M{s I H'{S) = Q}. 

The following quantity is often called the fractal dimension of the set S. It can 
be approximated in practical applications, by estimating the slope of the graph of 
the logarithm of the number of "boxes" of side length S that intersect S, versus the 
logarithm of S. 

Definition 4. The box-counting dimension of the set S* C X is defined to be 

dims J = hm - — 

s^o log(l/(5) 

if and only if this limit exists, where Ns{S) is the smallest number of sets of diameter 
(5 > that can cover S . 

In order to provide a precise calculation of box-counting and Hausdorff dimension 
of the attractor of an IFS we need the following condition. 

Definition 5. The IFS F is said to obey the open set condition if there exists a 
non-empty open set O such that 

M 

F{0) = U fm{0) C O, 
m— 1 

and 

fUO)nfi{O)^0 ifm^l. 

The following theorem provides the Hausdorff dimension of the attractor of an 
IFS in some special cases. 

Theorem 4. Let the IFS F consist of similitudes, that is fmix) = SmOmX + tm 
where Om is an orthonormal transformation on M.^ , Sm G (0,1), and tm & R^- 
Also let F obey the open set condition, and let A denote the set attractor of F. 
Then 

dim// A — dims A ~ D 
where D is the unique solution of 

M 

(2.6) ^™ - 1- 

m— 1 

Moreover, 

< H^{A) < oo. 

Proof. This theorem, in essence, was first proved by Moran in 1946, [221 ■ A full 
proof is given in ^21 P-118. | 

A good choice for the probabilities, which ensures that the points obtained from 
the random iteration algorithm are distributed uniformly around the set attractor 
in case the open set condition apphes, is pm = s^^. Note that the choice of D 
in Equation 12.61 is the unique value which makes {pi,p2, ...,pm) into a probability 
vector. 
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2.3. Code Space. A good way of looking at an IFS F as in H2.1|l is in terms of 
the associated code space E — {1, 2, M}°°. Members of S are infinite sequences 
from the alphabet {1, 2, M} and indexed by N. We equip S with the metric 
defined for u ^ xhy 



where k is the index of the first symbol at which uj and >f differ. Then (E, ds) is a 
compact metric space. 

Theorem 5. Let A denote the set attractor of the IFS F. Then there exists a 
continuous onto mapping _F : E — > A, defined for all (Ti(T20'3--- £ by 

F{ai<T2cr3...) = lim /^^ o /^^ o ... o fa^{x). 

k — >oo 

The limit is independent o/x G X and the convergence is uniform in x. 
Proof. This result is contained in yjj Theorem 3.1(3). | 

Definition 6. The point aia2<y3--- G is called an address of the point F {a ia2cr 3...) G 
A. 

Note that F : S ^ A is not in general one-to-one. 

The following theorem characterizes the measure attractor of the IFS F as the 
push-foward, under F : E ^ A, of an elementary measure p e P(E), the measure 
attractor of a fundamental IFS on E. 

Theorem 6. For each rn G {1, 2, M} define the shift operator s„j : E ^ E by 

s,n{o-ia2cr3...) = maia2cr3... 

\/aia2<T3... E E. Then s„i is a contraction mapping with contractivity factor . 
Consequently 

S -.^ {E; S1S2, sm;pi,P2, ■■■,Pm} 

is an IFS. Its set attractor is E. Its measure attractor is the unique measure n £ 
P(E) such that 

Tt{uJiUJ2UJ3-- e E|cJi = (Ti, CJ2 = 0-2, Wfe = CTfe} = • Prra ' ■•■ ' P<7k 

VfceN, V(7i,ai,...,(7fee{l,2,...,M}. 

If p, is the measure attractor of the IFS F, with F : A defined as in Theorem 
then 

fi^Fin). 

Proof This result is W Theorem 4.4(3) and (4). | 

We call S the shift IFS on code space. It has been well studied in the context 
of information theory and dynamical systems, see for example \7\, and results can 
often be lifted to the IFS F. For example, when the IFS is non-overlapping, the 
entropy (see |S] for the definition) of the stationary stochastic process associated 
with F is the same as that associated with the corresponding shift IFS, namely: 

-J^Pm^OgPrn. 
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3. Trees of Iterated Function Systems and Random Fractals 

3.1. SuperlFSs. Let {X,dx) be a compact metric space, and let M and N be 
positive integers. For n £ {1, 2, N} let F" denote the IFS 

F"={X;/r,/«,.../r,;pr,p^,...A} 

where each : X X is a. Lipshitz function with Lipschitz constant < Z < 1 
and the pj^ 's are probabihties with 

M 
m—l 

Let 

(3.1) T = {X; F\F\ F^ ; F,, P2, Pat}, 

where the P^s are probabihties with 

N 

Y,Pn = l,Pn>0 Vn 6 {1,2,..., AT}. 

n=l 

P„ > Vn e {1, 2, ...,N} and J2^^^ P„ = 1. 

As we will see in later sections, given any positive integer V we can use the set 
of IFSs J- to construct a single IFS acting on H(X)^. In such cases we call J- a 
superlFS. Optionally, we will drop the specific reference to the probabilities. 

3.2. Trees. We associate various collections of trees with !F and the parameters 
M and N. 

Let T denote the (M-fold) set of finite sequences from {1,2,...,M}, including 
the empty sequence 0. Then T is called a tree and the sequences are called the 
nodes of the tree. Foi i = iii2---ik G T let |i| — k. The number k is called the level 
of the node a. The bottom node is at level zero. If j = jij2---ji G T then ij is 
the concatenated sequence iii2---ikji32---ji- 

We define a level-k (M-fold) tree, or a tree of height k, to be the set of nodes 
of T of level less than or equal to k. 

A labelled tree is a function whose domain is a tree or a level-fc tree. A limb 
of a tree is either an ordered pair of nodes of the form {i,im) where i £ T and 
TO G {1,2,...,M}, or the pair of nodes (0,0), which is also called the trunk. In 
representations of labelled trees, as in Figures [TEl and [T^ limbs are represented by 
line segments and we attach the labels either to the nodes where line segments meet 
or to the line segments themselves, or possibly to both the nodes and limbs when 
a labelled tree is multivalued. For a two-valued labelled tree t we will write 

r(i) — {T{node i),T{limb i)) for i (zT, 

to denote the two components. 

A code tree is a labelled tree whose range is {1, 2, N}. We write 

n^{T I r:r-.{l,2,...,iV}} 

for the set of all infinite code trees. 

We define the subtree t : T {1,2, ...,N} of a labelled freer : T {1,2, ...,N}, 
corresponding to a node i — iii2---ik G T, by 

r(j) = T{ij) y j eT. 
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Figure 16. Pictorial representation of a level-4 2-fold tree labelled 
by the sequences corresponding to its nodes. The labels on the 
fourth level are shown for every second node. The line segments 
between the nodes and the line segment below the bottom node 
are referred to as limbs. The bottom limb is also called the trunk. 

In this case we say that r is a subtree of r at level k. (One can think of a subtree 
as a branch of a tree.) 

Suppose T and a are labelled trees with |r| < |cr|, where we allow \a\ = oo. We 
say that a extends r, and r is an initial segment of a, if a and r agree on their 
common domain, namely at nodes up to and including those at level |r|. We write 

T ^ a. 

If r is a level-fc code tree, the corresponding cylinder set is defined by 

H = [T]n := {(J eft : T -< a}. 
We define a metric on Q by, for uj ^ 

if k is the least integer such that wii) ^ >c{i) for some i ^ T with |i| = k. Then 
(ri, c?n) is a compact metric space. Furthermore, 

diam{Q) = 1, 

and 

(3.2) dzam([r]) = ^ 
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whenever r is a level-A: code tree. 

The probabihties P„ associated with T in Equation H3.1|l induce a natural prob- 
abihty distribution 

on il. It is defined on cyhnder sets [r] by 

(3.3) p([r])= n ^-W- 

i<l»l<kl 

That is, the random variables r(i), with nodal values in {1, 2, A''}, are chosen 
i.i.d. via Pr(r(j) — n) — Pn. Then p is extended in the usual way to the tr-algebra 
M{Q) generated by the cylinder sets. Thus we are able to speak of the set of random 
code trees f2 with probability distribution p, and of selecting trees cr G according 
to p. 

A construction tree for is a code tree wherein the symbols 1,2, and N are 
replaced by the respective IFSs F^,F^, and . A construction tree consists of 
nodes and limbs, where each node is labelled by one of the IFSs belonging to T. We 
will associate the M limbs that lie above and meet at a node with the constituent 
functions of the IPS of that node; taken in order. 

We use the notation JF(ri) for the set of construction trees for T. For cr e we 
write .F((t) to denote the corresponding construction tree. We will use the same 
notation !F{ct) to denote the random fractal set associated with the construction 
tree T{a), as described in the next section. 

3.3. Random Fractals. In this section we describe the canonical random fractal 
sets and measures associated with T in H3.1|l . 
Let !F be given as in (|3.1|l . let fc e N, and define 

Tk-.n^i M(X) ^ H(X), 

by 

(3.4) T,{a){K)^ U o /f ^) o ... o 

{ieT\ \i\=k} 

V (T G and K G ]HI(X). (The set Tk{<y){K) is obtained by taking the union of the 
compositions of the functions occurring on the branches of the construction tree 
J-{a) starting at the bottom and working up to the fc"^ level, acting upon K.) 

In a similar way, with measures in place of sets, and unions of sets replaced by 
sums of measures weighed by probabilities, we define 

J^fe : X P(X) ^ P(X) 

by 
(3.5) 

{i£T\ \i\ = k} 

V a G r2 and <^ G P(X) . Note that the Tk (c) (<r) all have unit mass because the 
sum (over m) to unity. 

Theorem 7. Let sequences of Junctions {Tk} and {Tk} be defined as above. Then 
both the limits 



T{a) = lim {Tk{(j){K)}, and T{a) = lim {Tk{cT){^)}, 

k — >-oo fc— *oo 
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exist, are independent of K and <;, and the convergence (in the Hausdorff and 
Monge Kantorovitch metrics, respectively,) is uniform in a, K, and<;. The resulting 
functions 

T : H(X) and f : P(X) 

are continuous. 

Proof. Make repeated use of the fact that, for fixed cr £ fJ, both mappings are 
compositions of contraction mappings of contractivity /, by Theorem 1. | 

Let 

(3.6) 9) = {J'ia) G H(X)|a e n}, and ^ = {T{cr) £ P(X)|ct G n}. 
Similarly let 

(3.7) «P = T{p) =po jP-i e P(^), and <p = T{p) = po T'^ e P(^). 

Definition 7. The sets ^ and are called the sets of fractal sets and fractal 
measures, respectively, associated with T . These random fractal sets and measures 
are said to be distributed according to and respectively. 

Random fractal sets and measures are hard to compute. There does not appear 
to be a general simple forwards (random iteration) algorithm for practical compu- 
tation of approximations to them in two-dimensions with affine maps, for example. 
The reason for this difficulty lies with the inconvenient manner in which the shift 
operator acts on trees cr e relative to the expressions (|3.4|l and (|3.5|) . 

Definition 8. Both the set of IFSs {F" : n = 1, 2, .., N} and the superlFS T are 
said to obey the (uniform) open set condition if there exists a non-empty open set 
O such that for each n G {1, 2, .., TV} 

M 

F"{0) = U C(0) C O, 

m—l 

and 

f"n{0) n fPiO) = 0W m,l G {1, 2, .., M} with m^l. 

For the rest of this section we restrict attention to (X, d) where X C and d is 
the Euclidean metric. The following theorem gives a specific value for the Hausdorff 
dimension for almost all of the random fractal sets in the case of " non-overlapping" 
similitudes, see Jl] and pT]. 

Theorem 8. Let the set of IFSs {F" : n — 1,2, ..,7V} consist of similitudes, 
i.e. f"i{x) — sl]^O^j^x + where 0"j is an orthonormal transformation on R^, 
s;^ G (0, 1), and C e M^, for all n G {1, 2, iV} and m G {1, 2, ...M}. Also let 
: n = 1,2, ..,N} obey the uniform open set condition. Then for *P -almost all 

dim// A — dims A — D 
where D is the unique solution of 

N M 
n— 1 7n—l 

Proof. This is an application of ^21 Theorem 15.2, p. 230. | 
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4. Contraction Mappings on Code Trees and the Space ftv 

4.1. Construction and Properties of ily- Let y e N. This parameter will 
describe the variability of the trees and fractals that we are going to introduce. Let 

^nxQx ...xn. 

" V ' 

V TIMES 

We refer to an element of as a grove. In this section we describe a certain IFS 
on ^}y , and discuss its set attractor fly- its points are (^-tuples of) code trees that 
we will call V -groves. We will find it convenient to label the trunk of each tree in 
a grove by its component index, from the set {1,2, V}. 

One reason that we are interested in Q.v is that, as we shall see later, the set of 
trees that occur in its components, called V -trees, provides the appropriate code 
space for a V-variable superfractal. 

Next we describe mappings from OX to that comprise the IFS. The mappings 
are denoted by 77" : Vl^ — » Vl^ for a e ^ where 

(4.1) ^:={{l,2,...,iV}x{l,2,...,n*'}''- 
A typical index a ^ A will be denoted 

(4.2) a = (ai, 02, .., ay) 
where 

(4.3) Ot, = (ni,;w„_i,t>„^2, •■•,Wu,m) 

where e {1,2, ...,7V} and e {1,2,...,]/} for m e {1,2,...,M}. 

Specifically, algebraically, the mapping 77° is defined in Equation H4.8|l below. But 
it is very useful to represent the indices and the mappings with trees. See Figure lTTI 
Each map 77" in Equation (|4.8|l and each index a € A may be represented by a V- 
tuple of labelled level-1 trees that we call (level-1 ) function trees. Each function tree 
has a trunk, a node, and M limbs. There is one function tree for each component of 
the mapping. Its trunk is labelled by the index of the component v S {1,2,..., V} to 
which it corresponds. The node of each function tree is labelled by the IFS number 
rti, (shown circled) of the corresponding component of the mapping. The m"' limb 
of the v*-^ tree is labelled by the number Vy^m G {1, 2, V}, for m G {1, 2, M}. 
We will use the same notation 77" to denote both a T^-tuple of function trees and 
the unique mapping 77" : — > Vl^ to which it bijectively corresponds. We will 
use the notation a to denote both a l^-tuple of function trees and the unique index 
a G ,4 to which it bijectively corresponds. 

Now we can describe the action of 77" on a grove uj g . We illustrate with 
an example where = 3, TV = 5, and M = 2. In Figure Hlsl an arbitrary grove 
UJ — (wi, 1^2,^3) e fl^ is represented by a triple of coloured tree pictures, one blue, 
one orange, and one magenta, with trunks labelled one, two, and three respectively. 
The top left of Figure ITSl shows the map 77", and the index a G A, where 

(4.4) 77^(^1, u;2,W3) = (^i(LJi,W2),C5(t^3,t^2),^4('^3,t^l)), 

and 

a = ((l;l,2),(5;3,2),(4;3,l)), 

represented by function trees. The functions {^„ : n — 1,2, ...,5} are defined in 
Equation l|4.7|l below. The result of the action of on uj is represented, in the 
bottom part of Figure El by a grove whose lowest nodes are labelled by the IFS 
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'1,1^ 




V 



Figure 17. Each map 77° in Equation (|4.8|) and each index a e A 
may be represented by a 1^-tuple of labelled level- 1 trees that we 
call (level-1 ) function trees. Each function tree has a trunk, a 
node, and M limbs. The function trees correspond to the compo- 
nents of the mapping. The trunk of each function tree is labelled 
by the index of the component v G {1,2,...,^} to which it cor- 
responds. The node of each function tree is labelled by the IPS 
number (shown circled) of the corresponding component of the 
mapping. The m*'' limb of the v^^ tree is labelled by the domain 
number w„.„j G {1, 2, V}, for m £ {1, 2, ...A/}. 



numbers 1, 5, and 4, respectively, and whose subtrees at level zero consist of trees 
from Lo located according to the limb labels on the function trees. (Limb labels of 
the top left expression, the function trees of 77°, are matched to trunk labels in the 
top right expression, the components of oj.) In general, the result of the action of 
77° in Figure [T7I on a grove lo E fl^ (represented by V trees with trunks labelled 
from 1 to N) is obtained by matching the limbs of the function trees to the trunks 
of the y-trees, in a similar manner. 

We are also going to need a set of probabilities {T'^ja £ A}, with 

(4.5) ^ T'" 1, T'" > Va e A 

These probabilities may be more or less complicated. Some of our results are 
specifically restricted to the case 

P P P 

(4.6) V r-„,r-„„...r-„ 



ni ^ n2 ■■ riv 



MV 



V 

which uses only the set of probabilities {Pi, P2, Pn} belonging to the superlFS 
. This case corresponds to labelling all nodes and limbs in a function tree inde- 
pendently with probabilities such that limbs are labeled according to the uniform 
distribution on {1,2,..., V}, and nodes are labelled j with probability Pj. 
or each n G {1,2,..., N} define the n*^ shift mapping ^„ : fi^^ — * J7 by 

(4.7) (e„(cc))(0) = n and (e„(cc;))(m^) - iu^ii) V i G T, to G {1, 2, M}, 

for u! — {llJi,uj2, ...,ujm) G f^*^- That is, the mapping ^„ creates a code tree with 
its bottom node labelled n attached directly to the M trees wi, W2, ■■■,(^m- 
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Figure 18. The top right portion of the image represents a grove 
uj ^ {uji,ll)2,uj3) £ il'^ hy Si triple of coloured tree pictures. The top 
left portion represents the map in Equation H4.4(l using level- 1 
function trees. The bottom portion represents the image r]°'{uj). 

Theorem 9. For each a E A define -q^ : Vl^ OX hy 



(4.8) 



'?n2('^i'2,iJ^' 



)i ^nv (^"v.i I ^fv,2 J ■■■7 '^f v,Af )) 



T/ien : fi^ — ^ 51^ is a contraction mapping with Lipshitz constant . 

Proof. Let a E A. Let {lji,u}2^ wy) and (cji.a;2...., cjy) be any pair of points in 
n^. Then 

d^v{r]°'{uji,uj2, ...,a;v'),?7°('^i,^2,---,wv)) = 
max do(C„„(wi,„ 1, 2 7 m)>C„J'^i>„ i^t^t, ...,a;„„ 

«G{l,2,...,V'} 

M v^{i.,2,...y} 

l\i i;e{l,2,...,y} M 



It follows that we can define an IFS $ of strictly contractive maps by 
(4.9) $ = {^7^;^7^P^a e y^}. 

Let the set attractor and the measure attractor of the IFS $ be denoted by fiy, 
and respectively. Clearly, Vly G ]HI(ri^) while G P(fi^). We call fiy the 
space of V-groves. The elements of Jly are certain T^-tuples of A/-fold code trees 
on an alphabet of N symbols, which we characterize in Theorem llll But we think 
of them as special groves of special trees, namely F-trees. 
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For all V G {1, 2, V^}, let us define Vlv,v C Vlv to be the set of w*'' components 
of groves in Q,v- Also let py G P(il) denote the marginal probability measure 
defined by 

(4.10) pv{B) := py{B, VL, n, n)VB e M{n). 

Theorem 10. For all v G {1, 2, .... V} we have 

(4.11) Qy^y = riyi {set of all V-trees). 

When the probabilities {P^ja G A} obey Equation J^.d] ), then starting at any initial 
grove, the random distribution of trees cj G O that occur in the v*^ components 
of groves produced by the random iteration algorithm corresponding to the IFS 
after n iteration steps, converge weakly to py independently of v, almost always, as 
n oo. 

Proof. Let S : f7 ^ — > ^ denote any map that permutes the coordinates. Then the 
IFS $ = {n^-^-q^^V^a G A} is invariant under S, that is $ = {D,^ ^Eri^E'^ , a G 
A}. It follows that EH.y — fly and S/iy = p-y It follows that Equation (|4.11|) 
holds, and also that, in the obvious notation, for any {Bi, B2, By) G (M>{Q))^ 
we have 

(4.12) pyiBi,B2,...,Bv)=PyiB^„Ba„...,B„,.). 
In particular 

Py{B) = py{B, n, n, ...,n) = py{fi, n, n, b, n, n) vs g M{n), 

where the "i?" on the right-hand-side is in the w*'* position. Theorem |51 tells us that 
we can apply the random iteration algorithm (Theorem[21) to the IFS $. This yields 
sequences of measures, denoted by {/iy : / = 1, 2, 3...}, that converge weakly to py 
almost always. In particular py\B,Q,,rt, ...,r2) converges to py almost always. | 

Let L denote a set of fixed length vectors of labelled trees. We will say that 
L and its elements have the property of V -variability, or that L and its elements 
are V -variable, if and only if, for all a; G L, the number of distinct subtrees of all 
components of uo at level k is at most V , for each level k of the trees. 

Theorem 11. Letuj G fi^. Thenui G fly if and only u! contains at mostV distinct 
subtrees at any level (i.e. lo is V-variable). Also, if a G il, then a is a V-tree if and 
only if it is V-variable. 

Proof. Let 

5 = {w G fl^ I [{subtrees of components of u) at level k}\ <V, V fc G N}. 

Then S is closed: Let {s„ G 5*} converge to s* . Suppose that s* ^ S. Then, at some 
level k € N, s* has more than V subtrees. There exists Z G N so that each distinct 
pair of these subtrees of s* first disagrees at some level less than I. Now choose n 
so large that s„ agrees with s* through level (fc -I- I) (i.e. dQv{sn,s*) < )■ 
Then s„ has more than V distinct subtrees at level k, a contradiction. So s* G S. 

Also S is non-empty: Let cr G 51 be defined by (j{i) = 1 for all i ^ T. Then 
(cr, (T, ...,cr) G S. 

Also S is invariant under the IFS $: Clearly any s G S* can be written s = ry°(s) 
for some a G A and sGS. Also, if s G S then 77°(s) e S. So S = U{r]"-{S) \ a G A}. 

Hence, by uniqueness, S must be the set attractor of $. That is, S — fly. This 
proves the first claim in the Theorem. 
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It now follows that if cr e il is a V-tree then it contains at most V distinct 
subtrees at level k, for each fc G N. Conversely, it also follows that if cr G has the 
latter property, then (cr, cr, cr) G fly, and so u G fiv,i- I 

Theorem 12. For all 

which implies that 

lim ^lv,i ~ ^, 

V—^oo 

where the convergence is in the Hausdorff metric. Let the probabilities {V^la G A} 
obey Equation J^.d] ). Then 

(4.13) c;r(n)(py,p) < 1.4 

which implies 

hm py = p, 

where p is the stationary measure on trees introduced in Section \3.^ and conver- 
gence is in the Monge Kantorovitch metric. 

Proof. To prove the first part, let M^^^ > V > . Let r be any level-fc code 
tree. Then r is clearly 1^-variable and it can be extended to an infinite l^-variable 
code tree. It follows that [r] n Vlv,i 7^ 0. The collection of such cylinder sets [r] 
forms a disjoint partition of by subsets of diameter jjk+r , see Equation (|3.2|l ). 
from which it follows that 

The first part of the theorem follows at once. 

For the proof of the second part, we refer to Section ^31 I 

Let Sy = A°°. This is simply the code space corresponding to the IFS <& defined 
in Equation (|4.9|l . From Theorem [3 there exists a continuous onto mapping $ : 
Sy ^ i7y defined by 

$(010203...) = lim rj"-^ o rj"-^ o ... o t]°'''(lu) 

k — ^oo 

for all 010203... G Sy, for any to. In the terminology of section [2.31 the sequence 
^10203.. G Sy is an address of the y-grove $(010203...) G f^y and Ey is the code 
space for the set of V-groves fiy. In general <& : Sy Oy is not one-to-one, as we 
will see in Section IT!^ 

4.2. Compositions of the Mappings 77°. Compositions of the mappings 77" : 
fl^ — > il^, o G ^, represented by V^-tuples of level-1 function trees, as in Figure 
1171 can be represented by higher level trees that we call level-k function trees. 

First we illustrate the ideas, then we formalize. In Figure El we illustate the 
idea of composing ^-tuples of function trees. In this example V — 3, N ~ 5, and 
M — 2. The top row shows the levcl-1 function trees corrresponding to o, 6, c G .4 
given by 

a = ((l;2,3),(3;l,3),(5;2,3)), 
6 =((4; 3,1), (2; 1,2), (3; 3, 2)), 
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Figure 19. Illustrations of compositions of function trees to pro- 
duce higher level function trees. Here V = 3, N = 5, and 
M ~ 2. We compose the level- 1 function trees corresponding to 
a = ((1; 2, 3), (3; 1, 3), (5; 2, 3)), b = ((4; 3, 1), (2; 1, 2), (3; 3, 2)), and 
c = ((1; 1, 2), (5; 1, 3), (4; 2, 3)). The top row shows the separate 
level-1 function trees, a, 6, and c. The second row shows the level-2 
function tree aob, and the function tree c. The last row shows the 
level-3 function tree a o b o c. 



and 

c=((l;l,2),(5;l,3),(4;2,3)). 

The first entry in the second row shows the 3-tuple of level-2 function trees aob. The 
bottom bracketed expression shows the 3-tuple of level-2 function trees ao6oc.Then 
in Figure 1^ we have represented 77" o 77'' o ri'^{uj) for cu — (wi, ix'2, W3) S fl^. The 
key ideas are (i) aobo c can be converted into a mapping 7^"°''°'= : fi^ _^ f^Sand (ii) 

The mapping 77"°''°'^ (cj) is defined to be the 3-tuple of code trees obtained by at- 
tatching the tree cuy to each of the top limbs of each level-3 function tree in aobo c 
with label v for all v e {1, 2, 3} then dropping all of the labels on the limbs. 

Next we formalize. Let k £ N. Define a level-k function tree to be a level-fc 
labelled tree with the nodes of the first k — 1 levels labelled from {1, 2, N} and 
limbs (i.e. second nodal values) of all k levels labelled from {1, 2, V}. We define 
a grove of level-fc function trees, say g, to be a V^-tuple of level-fc function trees, 
with trunks labelled according to the component number, and we define Gk to be 

the set of such g. Let G :— [J Gk- We will refer to a component of an element of 

fceN 

G simply as a function tree. For g d G we will write \g\ — fc, where fc e N is the 
unique number such that g £ Gk- 
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Figure 20. Illustration of the composition of the three mappings, 
rfo rf" o rf = jyaoSoc applied to a; = (tJi, a;2, t^^a) G ■ See also 
Figure El and the text. 



Then, for all g = (31,52, e G, for aU v e {1,2, V}, 

g^{node i) £ {1, 2, iV} \f i e T with |j| < \g\ - 1, 

and 

g-uihmb i) £ {1, 2, V} y i e T with |i| < \g\ . 

For all g,h ^ G we define the composition goh to he the grove of {\g\ + \h\)-\evel 
function trees given by the following expressions. 

(g o h)y{node i) = gy{node i) \/ i £T with \i\ < \g\ — 1; 
{g o h)y{limb i) — g^ilimh ij'i i with |i| < \g\ ; 
{g o h)y{node ij) = hg^^i^^f, ,){node j) \/ ij e T with \i\ = |.g| , |j| < \h\ - 1; 
{g o h)y{limb ij) = hg^(u.mb i){limb j) V ij e T with \i\ = \g\ , |j| < \h\ . 
For all 5 e G we define rjS : ^ fl^ by 

(rf {Lu))v{i) — gv{node i) \/ i E T with |i| < \g\ — 1, 

and 

W{^))v{ij) = ^g^(iimb i){node j) \/ ij e T with |i| = \g\ , |j| > 0. 

This is consistent with the definition of 77" : — * il^ for a G ^, as the following 
theorem shows. We will write 77^ to denote both the mapping 77^ : fly and 
the corresponding unique l/-tuple of level-fc code trees for all g d G. 

Theorem 13. For all g, h e Gk we have 
(4.14) 779°'' or/'. 

It follows that the operation o between ordered pairs of elements of G is associative. 
In particular, if (ai, 02, at) G A'^ then 



r;"! O 77"" O ... O 7/"'= = ^aioa20...0Qt, ^ 



(4.15) 

Proof. This is a straightforward exercise in substitutions and is omitted here. 
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We remark that Equations H4.15|l and H4.14|l allow us to work directly with 
function trees to construct, count, and track compositions of mappings 77° G A. 
The space G also provides a convenient setting for contrasting the forwards and 
backwards algorithms associated with the IFS <f>. For example, by composing func- 
tion trees in such as way as to build from the bottom up, which corresponds to 
a backwards algorithm, we find that we can construct a sequence of cylinder set 
approximations to the first component wi of w G fly without having to compute 
approximations to the other components. 

Let Gk{V) C Gk denote the set of elements of Gk that can be written in the 
form tti o 02 o ... o flfc for some a = (ai, 02, a^) G A'^. (We remark that Gk{V) is 
V^-variable by a similar argument to the proof of Theorem llll ) Then we are able to 
estimate the measures of the cylinder sets [(TyOio^ao fc)^^] by computing the prob- 
abilities of occurence of function trees g G GkiV) such that (77^)1 = (^010020. ..oafc-j^ 
built up starting from level-0 trees, with probabilities given by Equation H4.6|) . as 
we do in the Section ^31 

The labeling of limbs in the approximating grove of function trees of level-/c of 

<i>(aia2 ) defines the basic T^- Variable dependence structure of <i>(aia2 ). We 

call the code tree of limbs of a function tree the associated dependence tree. 

The grove of code trees for $(aia2 ) is by definition totally determined by the 

labels of the nodes. Nevertheless its grove of dependence trees contains all informa- 
tion concerning its V- Variable structure. The dependence tree is the characterizing 
skeleton of Variable fractals. 

4.3. A direct characterization of the measure py. Lct{F„}J^g = {-F", Fy}] 
be a sequence of random groves of levcl-1 function trees. Each random function 
tree can be expressed as = (7V„", L;'(1), L;;(M)) where the iV^'s and L^^'s 
corresponds to random labellings of nodes and limbs respectively. 

We assume that the function trees {F"}, are independent with Pr{N^ = j) = Pj 
and Pr{L'^{m) = k) ^ l/V for any w, fc G {1,..,^}, n G N, j G {l,...,iV}, and 
m G {1, ...,M}. 

First the family {L^, ...,Ly} generates a random dependence tree, K : T ^ 
in the following way. Let K{%) = 1. If K{ii, ...,in) = j, for some 
j = 1, then we define K{ii, ...,z„,i„+i) = LJ'(z„+i). 

Given the dependence tree, let /j — if Nl = 

The following theorem gives an alternative definition of py: 

Theorem 14. Py([T]) = Pr{Ii — r(i), V \i\ < |r|), where {Ii}i^T is defined as 
above, and t is a finite level code tree. 

Proof. It is a straightforward exercise to check that (FqoFio- • •oFi:_i)i is a level- fc 
function tree with nodes given by {Ii}\i\<k~i^ and limbs given by {i^(i)}|j|<j,. 
Thus 

Pr{I, = T{i), V i with |z| < fc - 1) 

(4.16) = PriiFo o Fi o • • ■ o Ffe_.i)i(node i) = r(i), V i with \i\ < k - 1). 
For 

a = (ai,a2, ...,ak) G A''. 

let 

-pa _ 'pai'pa-l -pOfc 



RANDOM IFSS 



31 



denote the probability of selection of 

f]^ — 1]°"^ o r/"^ o ... o rj"" = rj" 
By the invariance of /ly 

a.eA'' 

Now let r be a level-(A: — 1) code tree. Then 

aeA'' 

E 

{atE-A*^ I {?7^)i (node i) — r(i), V i} 

From this and H4.16|l it follows that Pv/([r]) = Pr{Ii = T{i), V i with \i\ < 
4.4. Proof of Theorem ini Equation (|4rT3]l . 



Proof. We say that a dependence tree is free up to level k, if at each level j, for 
1 < J < fcj the M-' limbs have distinct labels. If V is much bigger than M and 
k then it is clear that the probability of being free up to level k is close to unity. 
More precisely, if F is the event that dependence tree of (Fq o Fi o • • ■ o Fk~i)i is 
free and V > , then 

Py{[T])^M{[r],^.-.^))= E v-^^y{{v'')-\{[T],n,...,n)) 

-«^n O-7 n -n fi-^ 

i=i ^ ' i=\ ^ ' i=i 

^1--^ E*+ E ^ + -+ E 
\ i=i i=i i=i 

> 1 - ^ (Af2 + Af* + ... + A-/^'') 

^^2(fc+i) 2Af2fc 

> 1 z ^ r > 1 . 

2F(Af2 - 1) - W 
In the last steps we have assumed A/ > 2. 

Let S be the event that (7^ai°a20...oafc-j^ _ ^_ Then, itsm^ </ie independence of 
the random variables labelling the nodes of a free function tree, and using Equation 
(|3.3|l. we see that 

Pr(5|F)= n V<^^p{[T]). 

{ieT\ l<|i|<fc} 

Hence 

Pv{[r]) = PrlS*) = Pr(i^) Pr(S'|F) + Pr(F^) Pr(5|i^^) 

9 Af^*^ 

<Pr(5|F) + Pr(F^)<p([r]) + ^^. 

Similarly, 

p^([r]) > Pr(F) Pr(5|i^) 
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2k 



2M 
3V 



Hence 
(4.17) 



\Pvi[r])-p{[r])\^ 



In order to compute the Monge Kantorovitch distance d-p(Q,){py , p), suppose / : 
R is Lipshitz with iyip / < 1, i.e. — /(ro)] < da(w, to) V lj, G fi. Since 

diam{n) = 1 we subtract a constant from / and so can assume I/I < i without 
changing the value of / fdp — J fdpy . 

For each level-A: code tree t ^ Tk choose some uJr G [t] C fi. It then follows that 



fdp- / fdp 









reTk \ 



fdp- / fdp 



E / (/-/K))dp- E /((/ - /K))rfpy + E 



Py(M)) 



2k 



2 Af3fe 



3V 



since diam [t] < -jgr+r from Equation H3.2|l . |/(a)^)| < i, Lip / < 1, G H , and 
using Equation 14.171 Choose x so that ^ = M^^. This is the value of x which 



minimizes ^ 



2 _|_ Af 



3y 



M 
2V 



M 

Choose k so that A: < a; < fc + 1 . Then 

3 

= 2* 



37 V m" 



1 



3Af / 



7 f 

Hence Equation H4.13|l is true. | 

5. SUPERFRACTALS 

5.1. Contraction Mappings on and the Superfractal Set S}v.i- 

Definition 9. Let y G N, let A be the index set introduced in Equation \4-l^ , let 
T he given as in Equation L'i.l]) . and let probabilities {"P^la G A} be given as in 
Equation i4.5{ l. Define 

by 

M M M 

(5.1) r(x) === ( u f-HK.,,j, U /™ (^-^,™)' - U f^nHK^vj) 

m— 1 m— 1 m— 1 

yK ^ {Ki,K2,...,Kv) GH^,VaG^. Let 

(5.2) := {H^;r,7'",a G 
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Theorem 15. J-y is an IFS with contractivity factor I. 

Proof. We only need to prove that the mapping /" : is contractive 

with contractivity factor I, y a G A. Note that, ~ (Ki, K2, Km), L = 

{Li,L2.,...,Lm) em^', 

M M 
m—1 m—1 

< max{dH(/;;(i^™),/;^(L„))} 

m 
m 

= l-d^M{K,L). 
Hence, V(Xi, i^s, ■■■,Kv), {Li,L2, Ly) e H^, 

d^v (r (Xi, if2, ...,Kv),riLi,L2, ...,Lv)) 

M M 

= max{dH(U /m (^-.,™)' U /m (^-....))} 

in—l m—1 

< max{; • d^^^^((K^^ ^,Ky^.^,...,K^^,^,), 

(-^I'li.l 7 ^Vv,2 J ■•■7 ^Vv,M ))} 

< / ■ d^, ((ifl,if2, (^1,^2, V)). 

I 

The theory of IFS in Section 2.1 appUes to the IFS J-y- It possesses a unique 
set attractor i^y S E[(H^), and a unique measure attractor S P(IHI^). The ran- 
dom iteration algorithm corresponding to the IFS J-y may be used to approximate 
sequences of points {V- tuples of compact sets) in Sjy distributed according to the 
probability measure 

However, the individual components of these vectors in Sjy, certain special sub- 
sets of X, are the objects we are interested in. Accordingly, for all v G {1, 2, V}, 
let us define ijy.v C H to be the set of w**^ components of points in S)y. 

Theorem 16. For all v G {1, 2, V} we have 

^y,v = ^y,i- 

When the probabilities in the superlFS Ty are given by {4-6}) , then starting from 
any initial V -tuple of non-empty compact subsets o/X, the random distribution of 
the sets K gM that occur in the v*^ component of vectors produced by the random 
iteration algorithm after n initial steps converge weakly to the marginal probability 
measure 

^y^i{B) <Py (B, H, H, . .. , H)VB G B(H) , 
independently of v, almost always, as n 00. 

Proof. The direct way to prove this theorem is to parallel the proof of Theorem llOl 
using the maps {/" : ^ | a e ^} in place of the maps {77" : fi^ ^ | 
a G A}. 

However, an alternate proof follows with the aid of the map J- : Jl— > M(K) 
introduced in Theorem [7| We have put this alternate proof at the end of the proof 
of Theorem ini | 
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Definition 10. We callSjv,i a superfractal set. Points inSjv,i are called V-variable 
fractal sets. 

Example 4. See Figure [2] This example is similar to the one in Section 11.21 
It shows some of the images produced in a realization of random iteration of a 
N=V=2. Projective transformations are used in both IFSs, 

1.629a; + 0.135?/- 1.99 0.505a; + 1.935y - 0.216 
^ -0.780a; + 0.864y - 2.569 ' 0.780a; - 0.864y + 2.569^' 

1.616a; - 2.758y + 3.678 2.151a; + 0.567?/ + 2.020 
' ^ 1.664X - 0.944y + 3.883' 1.664a; - 0.944y + 3.883 

1.667a; + .098y- 2.005 0.563a; + 2.064y - 0.278 
^ -0.773a; + 0.790y - 2.575 ' 0.773a; - 0.790y + 2.575^' 

1.470X- 2.193?/ + 3.035 1.212a; + 0.686?/ + 2.059 
" ' 2.432a: - 0.581y + 2.872 ' 2.432a; - 0.581?/ + 2.872 

One of the goals of this example is to illustrate how closely similar images can 
be produced, with "random" variations, so the two IFSs are quite similar. Let us 
refer to images (or, more precisely, the sets of points that they represent) such as 
the ones at the bottom middle and at the bottom right, as "ti-trees". Then each 
transformation maps approximately the unit square O := {{x,y) \ < x < 1,0< 
y < 1}, in which each ti-tree lies, into itself. Both f2{x,y) and fi{x,y) map ti- 
trees to lower right branches of ti-trees. Both fl{x,y) and fi{x,y) map ti-trees 
to a ti-tree minus the lower right branch. The initial image for each component, 
or "screen", is illustrated at the top left. It corresponds to an array of pixels of 
dimensions 400 x 400, some of which are red, some green, and the rest white. Upon 
iteration, images of the red pixels and green pixels are combined as in Example 
n The number of iterations increases from left to right, and from top to bottom. 
The top middle image corresponds to the fifth iteration. Both the images at the 
bottom middle and bottom left correspond to more than thirty iterations, and are 
representive of typical images produced after more than thirty iterations. (We 
carried out more than fifty iterations.) They represent images selected from the 
superfractal Sj2,i according to the invariant measure ^2,1- Note that it is the 
support of the red and green pixels that corresponds to an element of 9)2,1- Note 
too the "texture effect" , similar to the one discussed in Example 13 

By Theorem El there is a continuous mapping Tv ■ '^v that assigns to 

each address in the code space Sy = A°°a. V-tuple of compact sets in S)v But 
this mapping is not helpful for characterizing Sjy because J-y ■ Sjy is not 

in general one-to-one, for the same reason that <& : Sy fiy is not one-to-one, as 
explained in Section IT!^ 

The following result is closer to the point. It tells us in particular that the set 
of V-trees is a useful code space for T^-variable fractals, because the action of the 
IFS $ on the space of T^-tuples of code trees is conjugate to the action of the IFS 
J-y acting on V-tuples of compact sets. (We are concerned here with the mappings 
that provide the correspondences between T^-groves, and y-trees, on the one hand, 
and points and probability distributions on i^y and Sjy.i, on the other.) 



superlFS with M= 
specifically 

fli.x,y) ^ 
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Figure 21. Sequence of images converging to 2-variable fractals, 
see Example^] Convergence to within the numerical resolution has 
occurred in the bottom left two images. Note the subtle but real 
differences between the silhouettes of these two sets. A variant 
of the "texture effect" can also be seen. The red points appear to 
dance forever on the green ti-trees, while the ti-trees dance forever 
on the superfractal. 

Theorem 17. Let a E A and rj"- : fi^ — > fi^ he defined as in Theorem\^ Let 
/° : be defined as in Theorem \15[ Let T : ]HI(X) be the mapping 

introduced in Theorem^ Define T : fJ^^ (]HI(X)) by 

for all {lui,uj2, ■■■,0Jv) ^ ■ Then 

(5.3) ^iv'^i^)) = n^i^)) aeA,ujen^. 
Also 

(5.4) T{Qv) = iOi^ and T{VLvs) = ^v,i, 

where Q.v,v denotes the set of v^^ components of members of fly ■ Similarly, when 
the probabilities in the IFS Ty of Equation \5. are given by Equation i4.6}) we 
have 

(5.5) J^ipy) = ^y, and T{py) = ^yt, 

where py is the marginal probability distribution given by Equation 

Proof. We begin by establishing the key Equation (|5.3|l . Note that from Theorem 
13 for anyis:GH(X), 

(5.6) T{uj) = {T{iui),T{u2),...,T{Lov)) 
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= ( lim {J^k{iOi){K)}, lim {Tk{oJ2){K)},..., lim {J^kMiK)}). 

k — >oo k — *oo k — *oo 

The first component here exemplifies the others; and using Equation H3.4|l it can 

be written 

(5.7) lim {Tk{uJi){K)} 

k — >oo 

= Jim { U o o ... o 

{*eT| \i\=k} 

Since the convergence is uniform and all of the functions involved are continuous, 
we can interchange the lim with function operation at will. Look at 

nn^)) = r( lim {M^i){K)}, Hm {Tk{u;2)iK)}, lim {J^kM{K)}) 

k — >oo k^oo k — >oc 

= lim {riTkiiOi)iK),Tkiu;2){K),...,Tkiiov){Km. 
By the definition in Theorem 1 151 Equation 15 .11 we have 

M M M 

m— 1 m— 1 m— 1 

By equation H3.4|l the first component here is 

M 



m— 1 

M 

™=1 {ieT| |i|=fc} 

= ^k+l{ini{^vt,^,^v^,2,--,^vt,M)){K), 

where we have used the definition in Equation (|4.7|l . Hence 

r{T{uj)) = hm {{r{Tu{uJi){K),Tk{u2){K),...,Tr,{uJv){K))} = 

k^oc 

k — >oc 

Comparing with Equations H5.6|l and H5.7|l . we find that the right hand side here 
converges to J-'{ri'^{u!)) as fc ^ oo. So Equation H5.3|l is true. 
Now consider the set We have 

aeA aeA aeA 

It follows by uniqueness that !F{nv) must be the set attractor of the IFS J-y- Hence 
J-{ilv) = Sjv which is the first statement in Equation (|5.4(l . Now 

Tiilv,!) = {T{uJi)\{uJi,U!2, ...,a;y) £ ily} 
= first component of {{!F{uJi), !F{lu2), ^i^^v)) \ (wi, 1^2, t^^v) G f^y } 
= first component of J-{Q,v) — first component of 9jv = ^v,i, 
which contains the second statement in Equation H5.4|l . 

In a similar manner we consider the push-forward under J- : — > of the 
invariant measure fiy of the IFS $ = {$7^; ry", P", a £ A}. lF{^y) is normalized, 
i.e. T{tiy) e P(H^), because T{tiy){M^) = Hy{T-^{W^)) = iJiyi^^). We now 
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show that J^{iiy ) is the measure attractor the IFS J-y- The measure attractor of 
the IFS $ obeys 

Applying J- to both sides (i.e. constructing the push-fowards) we obtain 

H^^v) - HY. ^^"(a^v)) - E ^"^W W)) = E ^"/"(-^(/v)), 

aG.A aeA aeA 

where in the last step we have used the key Equation l|5.3|l . So J-{p,y) is the measure 
attractor of the IFS J-y- Using uniqueness, we conclude J-{^y ) = *Py which is the 
first equation in Equation 15. 5() . Finally, observe that, for all B e 18(11), 

= T{fiy){B, H, H, H) = iiy{T^^{B, H, H, H)) 
= pLy{{T-'^{B),T-^{m),T-^{W),...,T-'^{m))) (using Equation OH ) 
= iJLy{{T-^{B),D.,n,...,n)) (since J^^ : ^2 ^ H) 
= py(.7^^^(B)) (by definition g^UIl) = ^{Pv){B). 
This contains the second equation in Equation (|5.5|l . 

In a similar way, we obtain the alternate proof of Theorem llfil Simply lift 
Theorem [ini to the domain of the IFS Tv using T -.Vl^ . | 

The code tree $(aia2a3...) is called a tree address of the V-variable fractal 
T{aia2a'i...). 

The mapping T : VLv,i — > together with Theorem 1111 provides a character- 
ization of V-variable fractals as follows. At any "magnification" , any V-variable 
fractal set is made of V "forms" or "shapes" : 

Theorem 18. Let M G 9jv.i be any V-variable fractal set. Let e > be given. Then 
M is a finite union of continuous transformations of at most V distinct compact 
subsets of"K, and the diameter of each of these transformed sets is at most e. 

Proof. Choose n so that Z" < e. Note that 

s^v=\Jn^v)= U r^or^o...or"(i5v'). 

aeA (ai,a2,...,a„)e-4 

Hence, since (M, M, M) e Sjv it follows that there exists (Ki, K2, Kv) G S)v 
such that 

M e first component of [j o o ... o {Ki, K2, Ky). 

{ai,a2,...,a„)eA 

Each set in the union has diameter at most e. | 

Example 5. This example is similar to Example 0] with M=N=V=2. The goal 
is to illustrate Theorem 1181 Figure |52 shows, from left to right, from top to bot- 
tom, a sequence of six successive images, illustrating successive 2-variable fractals, 
corresponding to a superlFS of two IFSs. Each IFS consists of two projective trans- 
formations, each mapping the unit square □ into itself. The images were obtained 
by running the random iteration algorithm, as described in Section f 1.21 The initial 
image on each screen was a blue convex region contained in a 400 x 400 array repre- 
senting □ , and the images shown correspond to one of the discretized screens after 
forty, forty-one, forty- two, forty-three, forty- four, and forty- five iterations. The key 
features of the transformations can be deduced from the images. (For example. 
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Figure 22. A sequence of 2- variable fractal sets (roughly accurate 
to viewing resolution), corresponding to M=N=2. Can you spot 
the projective transformations? See Example 01 Each of these six 
images images exhibits "2-variability" : at several scales, each looks 
as though it is the union of projective transformations of at most 
two distinct sets. 

one of the transformations of one of the IFSs, interpreted as a mapping from one 
screen to the next, maps the top middle image to the top two objects in the top 
left image.) Each of these images, at several scales, looks as though it is the union 
of projective transformations of at most two distinct sets. 

Theorem 19. The set of V-variable fractal sets associated with the superlFS Ty 
converges to the set of fractal sets associated with the superlFS T introduced in 
Section \cl.^ that is, in the metric o/H(]HI(X)), 

(5.8) Iimi3y,i=i3. 

V~*oo 

Moreover, if the probabilities {7^"|a G A} obey Equation i4.6j^ , then in the metric 
o/P(H(X)) 

lim = 

where *p is the stationary measure on random fractal sets associated with the su- 
perlFS T . 

Proof. We have, using the mapping T : il ^ M{X), 

lim ^)v,i = hm J^(f7ya) (by Theorem [T7|l 

V— *oo V—>-oc 

= T{ lim Qv,i) (since : ^ M(X) is continuous by Theoremd) 

V—*oo 

= (by Theoremini) 
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= io (by Equation ^^). 
Similarly, using the mapping J-' : Q ]P(X), have 

lim ^v.i — lim ^{Pv) (by Theorem 1 17|l 

V — >oo ' V — ^oo 

= T{ lim py) (since T : Vl ^ P(X) is continuous by Theorem O 

V^oo 

= F{p) (by Theorem HHl 
= *P (by Equation (|T7|) 1. 

I 

5.2. Contraction Mappings on and the Superfractal Measures Sjv,i- 
Recall that P = P(X). Let P^ = (P(X))^ In this section we follow the same lines 
as in Section I^TI constructing an IPS using the individual IFSs of the super IFS J^, 
except that here the underlying space is P^ instead of H^. 

Definition 11. Let 1/ £ N, let A be the index set introduced in Equation {.j.l^ , let 
J- he given as in Equation \S.l]) . and let probabilities {V°'\a G A\ be given as in 
Equation Define 

by 

MM M 

(5.9) r(A.) = (E^'mVrK,^), Ep™/™ (z^-..™)'---' Ep™ /™'(/^^v,.)) 

m— 1 m— 1 m— 1 

M = (Mi,M2> •■■-Mv) G P^- Let 

(5.10) J^-y := {P^;r,7'",a e ^}. 
Theorem 20. J-y is an IFS with contractivity factor I. 

Proof. We only need to prove that the mapping /" : P^ P^ is contractive 
with contractivity factor I, y a A. Note that, V/^ = (/ij^, /i2i Mm); V — 

((Pi,(P2,...,<^A/) eP*^ 

M M 

dvC^Plrflilirn), E P'^nfliVm)) 
m—1 m— 1 

M 



m—1 



< I ■ max{dp(/i„,,(p„J} 

m 

= I ■ dpM (/i, if ). 

Hence, V(/Zi,/X2, (<^i,</'2> e P^, 

dpv{f''{n^,p.^, ...,fiy),f"-{ip^,ip2, ...,(py)) 

M M 

= max{dp( E /m (m^.,„J. E P™ -^^ ("Pv^.J)} 

m—1 m—1 

< max{; • d^^, ((/i 

(V'«„,l'V'„„,2'--''^f.,M))} 

< ^C^,v((Ml,M2>■•■>MM):(<Pl>'^27■••>VM))• 
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The set attractor of the IFS Ty is f)v S ]HI(P ), a subset of P^, a set of ^-tuples 
of probabihty measures on X. As we wih see, each of these measures is supported 
on a 1/- variable fractal set belonging to the superfactal hv,\- The measure attractor 
of the IFS Tv is a probability measure e P(P^), namely a probability measure 
on a set of ^-tuples of normalized measures, each one a random fractal measure. 
The random iteration algorithm corresponding to the IFS Ty may be used to ap- 
proximate sequences of points in i3y , namely vectors of measures on X, distributed 
according to the probability measure '^y . 

As in Section [sm we define fiv,v to be the set of w*'' components of sets in ioy, 

for u e {i,2,...,y}. 

Theorem 21. For all v G {1, 2, V} we have 

When the probabilities in the IFS Ty are given by Equation i4.(^ , then starting at 
any initial V -tuple of probability measures on X, the probability measures fi G P(X) 
that occur in the v*^ component of points produced by the random iteration algorithm 
after n steps converge weakly to the marginal probability measure 

^y,i{B) ■■= ^v{B, P, P, P)VB G B(P), 
independently of v, almost always, as n oo. 

Proof. The direct way to prove this theorem is to parallel the proof of Theorem llOl 
using the maps {Z'' : P^ ^ P^ | a G A} in place of the maps {yy*^ : ^ fl^ \ 
a€A}. 

It is simpler however to lift Theorem EH to the domain of the IFS J-y using 
: P^ which is defined in Theorem [23 with the aid of the mapping 

!F : P = P(X) introduced in Theorem |7| We omit the details as they are 
straightforward. | 

We call Sjv,i a superfractal set of measures (of variability V). Points in Sjy^i are 
called V-variable fractal measures. 

Example 6. See Figure This example corresponds to the same superlFS as in 
Example^] The probabilities of the functions in the IFSs are Pi = Pi = 0.74, and 
pi ~ = 0.26. The IFSs are assigned probabilities Pi = P2 = 0.5. 

The following Theorem tells us in particular that the set of V^-trees is a useful 
code space for variable fractal measures, because the action of the IFS <& on the 
space of ^/-tuples of code trees is conjugate to the action of the IFS J-y acting on 
V-tuples of normalized measures. 

Theorem 22. Let a ^ A and r]" = Vt^ ^ be defined as in Theorem Let 
fa .pv ^pv defined as in Theorem\M LetT'.n^F^ P(X) be the mapping 
introduced in Theorem^ Define T : P^ = (P(X))^ by 

T{uji,uj2, ...,ujy) = {T{uJi),T(u}2), ...,T{ujv)), 
for all {uji,UJ2, ■■■jOJy) & ■ Then 



RANDOM IFSS 



41 




Figure 23. Three successive 2- variable fractal nieasures com- 
puted using the random iteration algorithm in Theorem l2 II applied 
to the superlFS in Example E| The pixels in the support of each 
measure are coloured either black or a shade of green, using a sim- 
ilar technique to the one used in Example [3 The intensity of the 
green of a pixel is a monotonic increasing function of the measure 
of the pixel. 



Also 

T{Qv) = and T{Qv,i) = ^v,i, 
where ^lv,v denotes the set of v*^ components of members of fly- Similarly, when 
the probabilities in the IFS Ty of Equation I5.1f)\) . are given by Equation i4.6\ ) we 
have 

T{Hv) = ^v, and T{py) ^v,i, 
where py is the marginal probability distribution given by Equation ^4.10\) . 

Proof. The proof is entirely parallel to the proof of Theorem 1171 using J- in place 
of J- and is omitted. | 

Definition 12. The code tree $(010203...) is called a tree address of the V-variable 
fractal measure ^y(ai0203...). 

The mapping J- : flv,i ^v.i together with Theorem 1111 allows us to charac- 
terize V-variable fractals as follows: 

Tileorem 23. Let ]1 G Sjv,i be any V-variable fractal measure. Let e > be 
given. Then Jl is a finite weighted superposition of continuous transformations of 
at most V distinct normalized measures supported on compact subsets o/X, and the 
diameter of the support of each of these transformed measures is at most e. 

Proof. Choose n so that Z" < e. Note that 

Hence, since (/i, /I, /i) G S)v it follows that there exists (tui, Ci72, roy) G 9]v 
such that 

/I G first component of o /°'^ o ... o (n7i, n72, ■ntjy). 

(ai,a2,...,a„)gX 
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Figure 24. Three successive 2-variable fractal measures, in shades 
of blue. Illustrates the "shapes" and "forms" theorem. See Exam- 
ple □ 

Inspection of Equation (|5.9() shows that each of the measures in the set of measures 
on the right-hand-side here is as stated in the theorem. | 

Example 7. See Figure |^ This corresponds to the same superlFS as used in 
Example 13 but here the measure is rendered in shades of blue to provide a pictorial 
illustration of Theorem|231 The three successive images were computed with the aid 
of the random iteration algorithm in Thcorcm l21l a new rendered measure theoretic 
image being produced at each iteration. At each discernable scale, approximately, 
each picture appears to have the property that it a superposition of a number of 
"little pictures" belonging to one of two equivalence classes. Pictures belonging to 
an equivalence class in this case are related by a projective transformation together 
with a scaling of brightness. 

Theorem 24. The set of V-variable fractal measures associated with the superlFS 
Tv converges to the set of fractal measures introduced in Section lS.Sl that is, in the 
metric o/H(P(X)) 

lim :?)vi ~ 

If the probabilities {P^ja G A} obey Equation i4.(^ , then in the metric o/P(P(X)) 

lim ^v,i - 

where *p is the stationary measure on fractal sets associated with the superlFS . 
Proof. We have, using the mapping T : ft ^ F = P(X), 

hm = lim ^(f^yj) (by Theorem !^ 

V—>oo y^oo 

T{ hm nv,i) (since T -M^F = P(X) is continuous by ThcorcmEl 

V—*OQ 

= (by Theorem 

^Sj (by Equation (E^J). 
We have, using the mapping T : ^ P(P), 

hm = lim $(Py) (by Theorem !^ 
~ !F{ lim py ) (since !F \ ^ ^ P(P) is continuous by Theorem [TJ 



RANDOM IFSS 



43 



= T[p) (by Theorem 1121) 
= $ (by Equation (j22|)). 

I 

5.3. Fractal Dimensions. Here we quantify and compare the Hausdorff dimen- 
sions of fractals corresponding to a (super) IFS of similitudes on for some 
if S N that obeys the open set condition in the following four cases: determin- 
istic fractals, standard random fractals, homogeneous random fractals {V = 1), 
and variable fractals {V > 1). The functions of the IFS are of the form 
fmi^) — ^mOm^ + tm where is an orthonormal transformation, £ (0, 1), 
and t^^ e M^, for all n e {1,2, ...,N} and m e {1,2, ...M}. 

5.3.1. Deterministic Fractals. In this case there is only one IFS, say F^. By Theo- 
rem 0] the Hausdorff dimension of the corresponding fractal set A is D, the unique 
solution of 

M 
m—1 

Example 8. Suppose K > 2. Let the IFS F^ consists of three similitudes with 
si = 32 = s^ — and that the fixed points are not coUinear. Then the set attractor 
of F^ is the Sierpinski triangle with vertices at the three fixed points. Its fractal 
dimension Di is given by 3 = 1 which implies Di = |^ = 1.585. 

Let the IFS consist of three similitudes with si — S2 — s^ = ^ and the same 
fixed points as F^ . Then the fractal dimension D2 of the set attractor of is the 
given by 3 = 1 which implies D2 — 1. 

5.3.2. Random Fractals. By Theorem [S] the Hausdorff dimension Dn of ?P-almost 
all of the random fractals sets for the superlFS T is given by 

N M 
n— 1 m—1 

Example 9. Let the superlFS be {D; F^ , F'^ ; Pi = P2 = 0.5} where the IFS's 
are defined in Example |H1 Then the fractal dimension Dj^ of *P-almost all of the 
random fractals in the set is given by i 3 + | 3 — 1 => — 1.262. 

5.3.3. Homogeneous Random Fractals (V = \). The case of homogeneous random 
fractals corresponds to 1^ = l.Each run of the experiment gives a different random 
Sierpinski triangle. 

Theorem 25. |15| . Let the superlFS T he as specified as in Theorem\^ Let V = \. 

Then for almost all A G fii.i 

AiTCLH A = D 

where D is the unique solution of 

N M 

Y,P„\nY,{s-,f ^l. 

n—l m—1 

Example 10. For the case of the superlFS in Example O whose 1-variable frac- 
tal sets we refer to as homogeneous random Sierpinski triangles, the Hausdorff 
dimension D of almost all of them is given by ilog(32TT) + llogl^^rj) = 0, 
=^d = 21og3/(log2-f log3) = 1.226. 
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5.3.4. V-Variable Fractals (V > 1). Let (ai,a2,...) G denote an i.i.d. se- 
quence of indices, with probabilities {■p"|a G A} given in terms of the probabihties 
{Pi, P2...., Py} according to Equation (|4.t)|) . Define, for a G [0, cx)) and a G A, the 
V X V flow matrix 

and let us write 

We think of sj^" as being the "flow" through the m*'* channel from screen v to 
screen w. where Vy„i = w. The sequence of random matrices M^^,... is 

i.i.d., again with probabilities induced from {Pi, P2,..., Py}. For any real square 
matrix M we define the norm ||Af|| to be the sum of the absolute values of its 
entries. By the Furstenberg Kesten Theorem \1'6\ . 

7(a) limk^^k-^ log ||Af ^(a) o • • ■ o M'=(a)|| 

exists and has the same value with probability one. Provided that the superlFS 
obeys the open set condition, we have shown in |H] that the unique value of D G 
[0, 00) such that 

l{D) = 

is the Hausdorff dimension of ^v.i almost all A G ^v,i- 

Kingman remarks that in general the calculation of 7 "has pride of place among 
the unsolved problems of subadditive ergodic theory", [211, P-897. However it 
is possible to estimate numerically. Namely, generate random copies of and 
iteratively compute M^, M^, ...,M^ and hence log ||M^(a) o • • • o M*'(a)|| for 
fc = 1, 2, ... The limit will give 7(0;). (Even for large V this will be quick since the 
are sparse.) One could now use the bisection method to estimate D. 

6. Applications 

Fractal geometry plays some role in many application areas, including the follow- 
ing. In biology: breast tissue patterns, structure and development of plants, blood 
vessel patterns, and morphology of fern fronds. In chemistry: pattern-forming al- 
loy solidification, and diffusion processes. In physics: transport in porous media, 
patterns formed during phase transitions in statistical mechanics, dynamical sys- 
tems, turbulence and wave propagation. In geology: particle size distribution in 
soil, and landscape habitat diversity. In computer science: digital image compres- 
sion and watermarking, compression of digital audio signals, image segmentation, 
and computer graphics. In engineering: wavelets, stochastic processes, rough sur- 
faces, antennae and frequency selective surfaces, stochastic optimal control, signal 
processing, and fragmentation of thin plates. 

In many of these areas it is clearly desireable to use random fractals; for example 
random fractals can be used in connection with diverse mathematical modeling ap- 
plication areas including Brownian motion, oil-wells, critical phenomena in statisti- 
cal physics, for example associated with lattice gasses and percolation, stock-market 
prices in finance, and in computer graphics they can be used to represent diverse 
picture types including natural images and textures. But random fractals are hard 
to compute, which may have held up the development of some applications, while 
deterministic fractals, which can be computed relatively easily, may not be rich 
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enough to provide convenient models for the applications to which one would want 
to apply them. 

Thus we believe that F-variable fractals could find many applications; they can 
be computed easily, with rapid access to many examples, contain a controllable 
amount of "randomness" , and have many of the advantages of fractals in general: for 
similitudes, with an open set condition, their fractal dimension may be computed, 
they are resolution independent, and they are in general geometrically complex 
at all levels of magnification, while being expressed with relatively small amounts 
of information, coefhcients of affine transformations and some probabilities, for 
example. 

6.1. Space-filling curves. Space-filling curves can be constructed with the aid of 
IFS theory, see for example 'T^, Chapter 9. These curves have many applications, 
including adaptive multigrid methods for numerical computation of solutions of 
PDEs and hierarchical watermarking of digital images. Here we note that inter- 
esting y-variable space-filling curves, and finite resolution approximants to them, 
can be produced. 

Example 11. Let M = 3,V = 2,N = 2. The IFS = {□; fl, f^Ji) consists of 
affine maps whose actions we explain with the aid of the left-hand diagram in Figure 
Esl □ is the unit square in the diagram, while /i (□) is the lower left square, /2 (D) is 
the upper left square, and /s (□) is the rectangle on the right. The transformations 
are chosen so that fl{OC) = OA, f^{OC) = AB, and f^{OC) = BC. Specifically 

The IFS = {□; , /|, /|} is explained with the aid of the right-hand diagram 
in Figure EHl fii^) is the lower left rectangle, /KD) is the upper left rectangle, 
and /|(n)_is the rectangle on the right; such that (OC) = OA', fi(OC) = WW, 
and /|(OC) = WC. Specifically fl{x,v) = (f y, ix), y) = + |, -ix -|- 

l)Ji{x,y) = {\x+l,-y+l). 

Neither of the IFSs here is strictly contractive, but each is contractive "on the 
average" , for any assignment of positive probabilities to the constituent functions. 
We assign probabilities Pi — P2 — 0.5 to the individual IFSs. An initial image 
consisting of the line segment OC is chosen on both screens, and the random 
iteration algorithm is applied; typical images produced after five iterations are 
illustrated in Figure !^ an image produced after seven iterations is shown in Figure 
123 Each of these images consists of line segments that have been assigned colours 
according to the address of the line segment, in such as way as to provide some 
consistency from one image to the next. 

6.2. Computer Graphics. New techniques in computer graphics are playing an 
increasingly important role in the digital content creation industry, as evidenced 
by the succession of successes of computer generated films, from "Toy Story" to 
"Finding Nemo". Part of the appeal of such films is the artistic quality of the 
graphics. Here we point out that T^-variable fractals are able to provide new types 
of rendered graphics, significantly extending standard IFS graphics 

Example 12. Here N=2, V=2, M=4. The two IFSs are given by 



{n;f-,f-,f-,f^;}n^{l,2}, 
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Figure 25. Transformations used for space-filling curves, see Ex- 
ample 




Figure 26. Low-order approximants to two 2- variable space filling 
curves, belonging to the same superfractal, see Example 1111 

where □ C M^, and each : □ ^ □ is a projective transformation. The colours 
were obtained as follows. A computer graphics rendering of the set attractor of 
is shown in Figure |2H1 and of in Figure 

The colouring of each of these two figures was obtained with the aid of an aux- 
iliary IFS acting on the cube C := [0,255]^ c given by G := {C ; gl' , , , g^} 
where each gm is a contractive (in the Euclidean metric) afhne transformation, 
represented by a 3 x 3 matrix and a 3 x 1 vector. For n G {1,2} discretized ap- 
proximations, of the same resolution, to the attractors of both IFSs F"- and Q were 
calculated via the deterministic algorithm (Corollary^; each pixel on the attractor 
of the IFS -F" was assigned the colour whose red, green, and blue components, each 
an integer from to 255, were the three coordinates of the point on the attractor 
of Q with the same code space address. At those points in the attractor of F" with 
multiple code space addresses, the lowest address was chosen. 

The superlFS we use is 

T = {□; F\F^; = 0.5, = 0.5} 

with V — 2. Then Figures l!!ni and mi show two examples, from among many different 
but similar ones, all equally visually complex, of computer graphics of 2-variable 
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Figure 27. Finite- resolution approximation to a 2-variable space- 
filling curve. See Example II II 

fractals for this superlFS, computed using the new random iteration algorithm. 
The images were rendered in much the same way as the images of the attractor 
sets of and were rendered above. The essential difference is the meaning of 
a "code space address" of a point on a l^-variable fractal, which we define to be the 
sequence of lower indices of a sequence of functions that converges to the point; for 
example, the point 

hm /I o /2 o /2 o /I o /2 o /I o ... o f^^^l (x) 

corresponds to the address 212122.. .mfc..., in the obvious notation. 

6.3. V-variable Fractal Interpolation. The technique of fractal interpolation 
has many applications including modelling of speech signals, altitude maps in geo- 
physics, and stock-market indices. A simple version of this technique is as follows. 
Let a set of real interpolation points in {{xi,yi) S = 0, 1, ...,/} be given. It 
is desired to find a continuous function / ; [a;o,a;/] R such that f{xi) — yi 
Vj e {0, 1, M}, such that its graph G = {{x,y) G : y = f{x)} is a fractal, 
possibly with specified fractal dimension. Introduce the IFS 

i^ = {K';/i,/2,...,/M} 

with 




Figure 29. The rendered set attractor of the IFS F'^m Example IT^ 




Figure 31. Another 2- variable fractal set for the superlFS IF in 
Example El 
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where the real coefficients am , e-m , Cm , d„i and Cm arc chosen so that 

frnixo, yo) = Vm-l, fmixo,yo) = Vm, 

and dm S [0, 1), for me {1,2, M}. Then the attractor of the IFS is the graph of 
a function / with the desired properties, its dimension being a function of the free 
parameters {dm : m = 1, 2, M}. 

Now let the superlFS ^ = {□;F\F2;pi = 0.5, = 0.5} for some V, consist 
of two IFSs both of which provide fractal interpolations of the data. Then all of the 
elements of the corresponding superfractal will be graphs of continuous functions 
that interpolate the data, have the property of ^-variability, and may be sampled 
using the random iteration algorithm. 

7. Generalizations 

It is natural to extend the notions of F-variable fractals, superlFS and super- 
fractal to include the case of maps contractive on the average, more than a finite 
number of maps, more than a finite number of IFSs, IFSs with a variable number 
of maps, IFSs operating on sets which are not necessarily induced by point maps, 
other methods of constructing the probabilities for a superlFS, probabilities that 
are dependent upon position etc. But for reasons of simplicity and in order to 
illustrate key features we have not treated these generalizations at any length. 
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